from scipy.special.orthogonal import p_roots

n=int(input("Input the order of Gauss quadrature:"))

[x,w]=p_roots(n)

print("Gauss points of {} order:".format(n))
print(x)
print("Gauss weigths of {} order:".format(n))
print(w)